Theory of the spatial structure of non-linear lasing modes 
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A self-consistent integral equation is formulated and solved iteratively which determines the 
steady-state lasing modes of open multi-mode lasers. These modes are naturally decomposed in 
terms of frequency dependent biorthogonal modes of a linear wave equation and not in terms of 
resonances of the cold cavity. A one-dimensional cavity laser is analyzed and the lasing mode is 
found to have non-trivial spatial structure even in the single-mode limit. In the multi-mode regime 
spatial hole-burning and mode competition is treated exactly. The formalism generalizes to complex, 
chaotic and random laser media. 



The steady-state electric field within and outside of 
a single or multi-mode laser arises as a solution of the 
non-linear coupled matter-field equations, the simplest of 
which are the two-level Maxwell-Bloch equations treated 
below. While the basic equations involved have been 
known for many years, and many aspects of their tem- 
poral dynamics have been studied relatively little 
progress has been made in understanding the spatial 
structure of the non-linear electric field, particularly in 
the case of multi-mode solutions for which spatial hole- 
burning and other non-linear effects are critical. It is 
natural to attempt to understand the non-linear solu- 
tions in terms of solutions of a linear wave equation. The 
two standard choices are either the hcrmitian solutions 
of a perfectly reflecting (closed) passive laser cavity Q, 
or the non-hermitian non-orthogonal resonances of the 
open passive cavity 0,0). In fact the intuitive picture 
of a lasing mode is that it arises when one of the res- 
onances of the passive cavity is "pulled" up to the real 
axis by adding gain to the resonator. Often compari- 
son of the numerically generated lasing modes with cal- 
culated linear resonances do show strong similarities in 
spatial structure, providing useful interpretation of lasing 
modes 0, , although not a predictive theory. However 
with the current interest in complex laser cavities based 
on wave-chaotic shapes El 01 , photonic bandgap media 
P. Hfl] or random media JlJ1ld| it is important to have 
a quantitative and predictive theory of the lasing states, 
as the numerical simulations required to solve the time- 
dependent Maxwell-Bloch equations are time-consuming 
and not easy to interpret. 

In recent work we have formulated a theory of steady- 
state multi-mode lasing which addresses these concerns 
^3 ■ The theory implies that the natural linear basis for 
decomposing lasing solutions is the dual set of biorthogo- 
nal states corresponding to constant outgoing and incom- 
ing Poynting vector at infinity at the lasing frequencies 
(referred to as "constant flux" (CF) states). Our the- 
ory shows that even in conventional lasers it is incorrect 



to regard the lasing modes as corresponding to a single 
resonance of the passive cavity and that multiple spa- 
tial frequencies occur even when there is a single lasing 
frequency close to the frequency of a single passive cav- 
ity resonance. These multiple spatial frequencies arise 
because several CF states contribute to a single lasing 
mode. Note that biorthogonal modes have been used ex- 
tensively in resonator theory (notably for the case 
of unstable resonators) but have not previously been ap- 
plied to multimode lasing theory. For multimode lasing 
the main difficulty is treating modal interactions and the 
related effects of spatial hole-burning 0| ■ We sketch be- 
low an efficient method for treating these effects exactly 
which can in principle be used in designing laser cavities 
to predict power output and tailor the mode spectrum of 
the laser. The techniques are illustrated for the simple 
case of a one-dimensional edge-emitting laser. 

We begin with the semiclassical laser equations within 
the rotating wave and slowly-varying envelope approxi- 
mations (see Ref. 01 f° r a derivation) 
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describing a laser comprised of a uniform gain medium 
of two- level atoms (with level spacing huj a ) embedded 
in a background dielectric medium/cavity with arbitrary 
spatially varying index of refraction n(x). Here e(x,t), 
p(x,t) are the envelopes of the field and polarization, 
D[x, t) is the inversion, Dq is the pump strength, g is the 
dipole matrix element, jj_ and 7|| are damping constants 
for p, D. As usual the fast variation of the fields at uj a 
is removed and the actual fields are given by (E, P) = 
(e,p) e l " at + c.c. For simplicity we take e, p to be scalar 
fields, appropriate for the ID case with TM polarization 
that we will discuss below; the same scalar form would 
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apply for planar random or chaotic cavities. n(x) is the 
(possibly) spatially dependent index of refraction of the 
cavity. 

We assume a steady-state lasing solution which is 
multi-periodic in time: e(x,t) = ^ M ^ni x ) e ~' if2 ' J *, 
p(x, t) = ^2 /1 P f j,(x)e~' infit . In contrast to standard modal 
expansions, not only the lasing frequencies but the 

spatial mode functions {^^(x)} are assumed to be un- 
known. 

Such a multi-periodic solution requires that the inver- 
sion is approximately stationary |l7|. implying [l3| that 
each lasing mode must satisfy the self-consistent equa- 
tion. 
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where g(^l^) is the gain profile evaluated at the las- 
ing frequency, D is the cavity domain and G is the 
Green function of the cavity wave equation [V 2 + 
n 2 (x)k 2 ] G(x, x'\oj) = 8 3 (x — x') with purely outgoing 
boundary conditions and k = lu/c is an external wavevec- 
tor of the lasing solution at infinity (for multimodc solu- 



tions k 



f2 M /c). Henceforth we set c = 1 and use 



wavevector and frequency interchangeably. With these 
non-hermitian boundary conditions the spectral repre- 
sentation G(x,x'\k) is of the form: 



G(x,x'\k)=J2 
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Here the functions (p m (x,k) are the CF states which 
satisfy — V 2 (p m (x, k) = n 2 (x)k 2 n ip m (x, k) with the non- 
hermitian boundary condition of only outgoing waves 
of wavevector k at the cavity boundary. For the spe- 
cial case of a ID cavity of length a considered below 
(see Fig. |2 inset) this condition is just d x ip m (x)\ a = 
+ikip rn (a). Note this differs subtly but importantly 
from the quasi-bound state boundary condition where 
the complex eigenvalue k m replaces the real wavevec- 
tor k The dual set of functions (p m {x,k) sat- 
isfy the complex conjugate differential equation with 
the boundary conditions d x (p m {x)\ a = —ik(p m (a) cor- 
responding to constant incoming flux. In general 
these functions satisfy the biorthogonality relations [lgj : 
J v dx <f^ n (x,k)n 2 (x)ip n (x,k) = rj m (k)5 mn , and are also 
complete. These relations make it possible to expand an 
arbitrary lasing solution 
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so that each vE^ is a vector in the space of biorthogonal 
functions. Here, Lp^x) = (p m (x,k fJi ) and in what fol- 
lows we define kl^ = k^k^). By substitution of © into 



Eq. |0J and use of the biorthogonality relations one finds: 
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where we have rescaled the pump 2irk a g 2 Do/h r y± — > Dq, 
and measured electric field in units e c = HyFffY±/2g. 

Equation (0 is the key result of our work; it deter- 
mines the lasing mode(s), each of which is a superposition 
of CF states which depends on its lasing frequency, fc^, 
and the pump power, Dq. It is useful to regard Eq. (jTJ 
as defining a map of the complex vector space of coeffi- 
cients a M = (a^ , Q-2 , . • .) into itself, where the actual non- 
zero lasing solution is a fixed point of this map. Above 
the lasing threshold for each mode, D^ t , we find that the 
non-zero solutions are stable fixed points and trial vectors 
flow to them under iteration of the map while the trivial 
zero solutions, which below D$ t arc stable, become un- 
stable. Note that the map is proportional to [k^ — k m ] , 
favoring the CF state with complex wavevector close to 
the real lasing wavevector, and it is also proportional to 
[7_L — ifcp] -1 , insuring that the lasing frequency is near 
the center of the gain profile. It can be shown that in the 
high-finesse limit in which the imaginary part of the CF 
frequency is very much smaller than the real spacing be- 
tween them only one CF state dominates the lasing state 
(the "single-pole approximation"), and this CF state is 
virtually identical to the corresponding linear resonance 
In this limit the picture of a single resonance being 
"pulled" up to the real axis is valid. However in many 
realistic cases this limit is not realized and the actual las- 
ing solution is the superposition of CF states determined 
byEq. J7J). 

Eq. 10 determines the lasing frequencies as well. For 
the first lasing mode and a uniform index resonator this is 
particularly simple. At threshold the CF states in the de- 
nominator of the integrand can be ignored and biorthog- 
onality leads to a simple relation: 
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For a non-trivial solution we must have a m =/= for some 
m and hence the coefficient must be real and equal to 
unity. The reality condition determines that the possible 

lasing frequencies at threshold are kjf m ^ = ' y± qm( ~ kt (m ] 

where k m = q m (k) — in m (k) (we suppress the index /x 
here). Furthermore, the modulus unity condition deter- 
mines the threshold pumping: 
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The CF state m and associated frequency leading to the 
lowest threshold will be the first lasing mode. Note that 
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in contrast to the traditional mode-pulling formula 0] 
where the cavity mode frequency is a fixed value, here the 
single-mode laser frequency is determined by the solution 
of a self-consistent equation (a transcendental equation 
for the ID case). Nonetheless for high finesse cavities 
this condition agrees with standard results: the lasing 
frequency is very close to the cavity resonance nearest 
to the gain center, pulled towards the gain center by an 
amount which depends on the relative magnitude of 7j_ 
vs. K m (which is approximately the resonance linewidth) 

Above threshold the lasing mode is found by initially 
choosing the lasing wavevector, k = k^ n \ calculating the 
sets {</? m },{<^ m } corresponding to that choice and then 
iterating Eq. JJJ) starting from a trial vector a(0) to 
yield output a(l). A natural choice for the initial vector 
is a m = 1, a p = 0, V p 7^ m, where m is the dominant com- 
ponent at threshold, calculated from the above relations. 
As noted, the "lasing map" has the property that below 
the lasing threshold, D t, the iterated vector, a — > 0, 
and above this threshold it converges to a finite value 
which defines the spatial structure of the lasing mode 
in terms of the CF states. There is one crucial addi- 
tion necessary to complete the algorithm. Note that Eq. 
(|7|) is invariant under multiplication of the vector a by a 
global phase e 10 , so iteration of (Q) can never determine a 
unique non-zero solution. Therefore it is necessary to fix 
the "gauge" of the solution by demanding that we solve 
Q with the constraint of a certain global phase (typi- 
cally we take the dominant a m to be real). Thus after 
each iteration of 10 we must adjust the lasing frequency 
to restore the phase of a m ; it is just this gauge fixing 
requirement which causes the lasing frequency to flow 
from our initial guess to the correct value above thresh- 
old. The invariancc of Q under global phase changes 
guarantees that the frequency thus found is independent 
of the particular gauge choice. For multi-mode lasing we 
repeat this procedure for each vector a M . In the single- 
mode regime only one of these vectors will flow to a non- 
zero fixed point. This behavior of the multi-mode lasing 
map is illustrated in Fig. ^ for the simple uniform index 
ID cavity corresponding to an edge-emitting laser with 
a perfect mirror at the origin and an index step at x = a 
(inset to Fig. |2J. Below the first threshold, determined 
from Eq. [3] above, the entire set of vectors a M flow to 
zero; above that threshold the first lasing mode turns on 
and its intensity grows linearly with pump strength. Due 
to its non-linear interaction with other modes, the turn- 
on of the second and third lasing modes is dramatically 
suppressed, leading to a factor of four increase in the 
interval of single-mode operation. The intensity shows 
slope discontinuities at higher thresholds as seen in nor- 
mal laser operation. Note that in this approach effects 
of spatial hole-burning and mode competition are treated 
exactly and not in the near-threshold a ppr oximation (cu- 
bic non-linearity) traditionally used which greatly 
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FIG. 1: Convergence and solution of the multimode lasing 
map for ID resonator with 710 = 1.5, k a — 19.89 and j± = 4.0 
vs. pump Do. Three modes lase in this range. At threshold 
they correspond to CF states m = 8, 9, 10 with threshold 
lasing frequencies fcj = 18.08, fc' 9 - 1 = 19.91, and fc^ 10 ' = 
21.76, and non-interacting thresholds D 9 ^ — 1.204, Dq 10 ^ = 
1.445, D^ = 1.482 (green dots). kf ] « k a and m = 9 has 
the lowest threshold. Due to mode competition, modes 2,3 
do not lase until much higher values (Do = 2.25, 2.53). Each 
mode is represented by an 11 component vector of CF states; 
we plot the sum of |a m | 2 vs. pump Do. Below threshold the 
vectors flows to zero (blue dots). For Do > 1.204 the sum 
flows (red dots) to a non-zero value (black dashed line), and 
above Do = 2.25, 2.53, two additional non-zero vector fixed 
points (modes) are found (convergence only shown for modes 
1,2). 

underestimates the output power [l3| . 

We now consider the spatial structure of the CF states 
defining the lasing modes. The linear resonances (quasi- 
bound states) of this system are easily found 01 1 they 
are complex sine waves of wavevector nok^a = n(m + 
1/2) - i/21n[(no + l)/(n - 1)]. The constant linewidth 
Im[km b ^] = -Km'' follows from the Fresnel transmissiv- 
ity of dielectric boundary at normal incidence. Note that 
Km'' > here, corresponding to amplification with in- 
creasing x. From Eq. JSJ we expect the lasing mode to 
involve several CF states and differ most from a single 
cavity resonance for a low finesse cavity; thus we con- 
sider relatively small index, no = 1.5. The CF states 
depend on the lasing wavevector, k. In Fig. (0) we 
choose it to be the real part of the wavevector of the 9 th 
resonance, kg b , and plot the 7 closest CF eigenvalues, k m 
(inset, Fig. [2). The 9 th CF state has k 9 w kf and within 
the cavity is very close to the m = 9 resonance [lflj . 
but the other k m , while they have Re[k rn ] sa Reftm^] 
(hence similar FSR) have substantially larger or smaller 
K m than Km''. Hence only the m = 9 CF state is close 
to a linear resonance, emphasizing that CF states are 
not resonances. We plot several of these states in Fig. |2 
showing their different amplification rates. The actual 
lasing mode will be the sum of several of these CF modes 
with different spatial frequencies and amplification rates. 
In Fig. we plot such a mode. Standard modal expan- 
sions in laser theory are equivalent to choosing only the 
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FIG. 2: The spatial structure of three consecutive CF modes 
calculated at the external frequency k = 19.90 and no — 
1.5. The CF frequencies are k m = 19.90 - 0.536 i (black), 
21.99 - 0.466 i (red), and 17.81 - 0.640 i (blue). The upper 
inset compares the CF mode frequencies (+) with the quasi- 
bound mode frequencies (o) in the range m = [6, 12]. Lower 
inset: schematic of the edge-emitting laser cavity. 

central CF state and missing the contribution of these 
spatial "side-bands" 01 . The inset to Fig. [3] shows that 
near threshold only one CF state dominates (one can 
show that the other components are of order the cube 
of the dominant component). But well above threshold 
the two nearest neighbor CF states are each 15% of the 
main component and since one of these has higher am- 
plification rate, the final effect is to increase the output 
power by more than 43% (see Fig. 0). The sidebands 
are still 6% of the dominant component when the index 
is increased to no = 3, leading to an increase in out- 
put power by 26% (see inset, Fig. The formalism 
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FIG. 3: Non-linear electric held intensity for a single-mode 
edge-emitting laser with no = 1.5, k a = 19.89, j± = 0.5, Do = 
9. The full field (red line) has an appreciably larger amplitude 
at the output x = a than the "single-pole" approximation 
(blue) which neglects the sideband CF components. Inset: 
The ratio of the two largest CF sideband components to that 
of the central pole for no — 1.5 (□, x) and no = 3 (□,+) vs. 
pump strength Do- 

we have presented here is ideal for treating random or 



wave-chaotic lasers for which the output directionality, 
output power and mode spectrum are very hard to pre- 
dict with heuristic arguments. Moreover in these systems 
the finesse typically is parametrically smaller than unity, 
suggesting that each lasing mode will consist of many 
CF states without a dominant component. We have con- 
sidered here the borderline case of a cavity with finesse 
of order unity, so there is a dominant CF component, 
but the spatial sidebands are still appreciable. Since our 
method treats the non-linearity and mode-competition 
exactly we anticipate that it may be useful in designing 
efficient semiconductor microcavity lasers. 
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